Confined Coulomb systems with adsorbing boundaries: 
the two-dimensional two-component plasma 
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Using a solvable model, the two-dimensional two-component plasma, we study a Coulomb gas 
confined in a disk and in an annulus with boundaries that can adsorb some of the negative particles 
of the system. We obtain explicit analytic expressions for the grand potential, the pressure and 
the density profiles of the system. By studying the behavior of the disjoining pressure we find that 
without the adsorbing boundaries the system is naturally unstable, while with attractive boundaries 
the system is stable because of a positive contribution from the surface tension to the disjoining 
pressure. The results for the density profiles show the formation of a positive layer near the boundary 
that screens the adsorbed negative particles, a typical behavior in charged systems. We also compute 
the adsorbed charge on the boundary and show that it satisfies a certain number of relations, in 
particular an electro-neutrality sum rule. 

Keywords: Coulomb systems, two-component plasma, adsorbing boundaries, soap films and bubbles, micelles 
and vesicles, disjoining pressure, charge density 
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I. INTRODUCTION 

In this paper we study the classical (i.e. non- 
quantum) equilibrium statistical mechanic properties of 
confined Coulomb systems with adsorbing boundaries. 
A Coulomb system is a system of charged particles inter- 
acting through the Coulomb potential. There are several 
interesting realizations of Coulomb systems with several 
applications such as plasmas, electrolytes, colloidal sus- 
pensions, etc... In the present paper we are interested 
in the case where the Coulomb system is confined with 
boundaries that can attract and adsorb some particles of 
the system. Our study of this kind of systems will be done 
using a solvable model of Coulomb system: the symmet- 
ric two-dimensional two-component plasma, a system of 
two kind of oppositely charged particles ±q at thermal 
equilibrium at an inverse temperature j3 = (ksT)^ 1 . The 
classical equilibrium statistical mechanics of the system 
can be exactly solved when Pq 2 = 2. 

One can think of several examples where the present 
situation of Coulomb systems confined with adsorbing 
boundaries is relevant, for instance in a plasma or an elec- 
trolyte near an electrode with adsorbing sites [l], |2| ■ An- 
other situation in which we will focus our attention is in 
solutions of amphiphathic molecules and ions for example 
in soap films and bubbles. Amphiphathic molecules have 
an hydrophobic tail and an hydrophilic charged head 
(usually negative) and, for this reason, when they are 
submerged in water they rearrange themselves in such a 
way as to minimize the contact of the hydrophobic tails 
with the surrounding medium. They can achieve con- 
figurations such as bilayers, micelles and vesicles, among 
others. 

In a previous paper Q we studied a soap film by model- 
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ing it as a Coulomb system confined in a slab. A soap film 
can be seen as a system of amphiphathic molecules in a 
bilayer configuration with a water inner layer. The over- 
all neutral system with negatively charged amphiphathic 
anions and positive micro-cations (usually Na + ) in wa- 
ter was modeled as a two-dimensional two-component 
plasma. The two dimensions were in the breadth of the 
film not on the surface: we studied a cross section of the 
film. Because of the hydrophobicity of their tails, the 
soap anions prefer to be in the boundaries of the film. 
This was modeled by a one-body attractive short-range 
external potential acting over them. This means that 
the negative particles felt an attractive potential over a 
small distance near each boundary. In this sense the neg- 
ative particles of the system can be "adsorbed" by the 
boundary. 

In Ref. Q we found exact expressions for the density, 
correlations and pressure inside the film. By studying 
the disjoining pressure, we were able to conclude that 
the Coulomb interaction plays an important role in the 
collapse of a thick soap film to a much thinner film. Ac- 
tually if a large number of amphiphathic molecules are in 
the boundary due to a large strength of the attractive ex- 
ternal potential near the boundary, the system is always 
stable. On the other hand if the attractive potential near 
the boundary is not strong enough a thick film will not 
be stable and will collapse to a thin film. 

These results can be compared (qualitatively due to 
the simplicity of the model under consideration) to the 
experimental situation of the transition of a thick film 
to a thin black film. These black film phenomena occur 
when the soap films width is smaller than visible light 
wavelength and it is seen black. Two types of black films 
are observed experimentally: the common black film and 
the much thinner Newton black film. 

It is interesting to know to what extent the results 
of our previous work |3( depend on the geometry. Here 
we study this two-component plasma system confined 
in two other geometries, a disk and an annulus. As 
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we mentioned before, amphiphathic molecules in water 
can achieve micelles and vesicles among others config- 
urations. In two dimensions, a cylindrical micelle can 
be seen as a disk and a cylindrical vesicle as an annu- 
lus. If the length of the cylindrical micelle or vesicle if 
much larger that its radius it is reasonable to assume 
that the system in invariant in the longitudinal direc- 
tion and so we study only a cross section of the sys- 
tem: a disk or an annulus. The Coulomb interaction 
is then the two-dimensional Coulomb potential which is 
v c (r) = — ln(r/d) for two particles at a distance r of each 
other. The length d is an arbitrary length which fixes the 
zero of the potential. 

Two-dimensional Coulomb systems with log interac- 
tion have properties that are similar to those in three 
dimensional charged systems with the usual 1/r poten- 
tial. They satisfy Gauss law and Poisson equation in two 
dimensions. Several universal properties, such as screen- 
ing effects, are direct consequences of the harmonic na- 
ture of the — ln(r/d) and 1/r potentials, which are the 
solutions of the two- and three-dimensional Poisson equa- 
tion. Therefore the exact solutions obtained for the 2D 
models play an important role in understanding real 3D 
Coulomb systems. 

The rest of this work is organized as follows. In sec- 
tion [n] we present in detail the system under consider- 
ation and briefly review how this model can be exactly 
solved. In section ITTT1 we compute the grand potential of 
the system and the disjoining pressure and study the sta- 
bility of the system. In section llVl we compute the density 
profiles of the different types of particles in the system 
and the adsorbed charge on the boundaries. Finally, we 
conclude recalling the main results of the present work. 



II. THE MODEL AND METHOD OF SOLUTION 

The system under consideration is a two-dimensional 
system composed of two types of point particles with 
charges ±q. Two particles with charges sq and s'q at 
a distance r apart interact with the two-dimensional 
Coulomb potential — ss'q 2 ln(r/d) where d is an arbi- 
trary length. This system is known as the symmetric 
two-dimensional two-component plasma. 

When the Boltzmann factor for the Coulomb potential 
is written, the adimensional coulombic coupling constant 
appears: T — q 2 /ksT = (3q 2 . Notice that in two di- 
mensions q 2 has dimensions of energy. For a system of 
point particles if Y > 2 the system is unstable against the 
collapse of particles of opposite sign, the thermodynam- 
ics of the system are not well defined unless one consid- 
ers hard-core particles or another regularization proce- 
dure (for instance a lattice model instead of a continuous 
gas). On the other hand, if Y < 2 the thermal agitation 
is enough to avoid the collapse and the system of point 
particles is well defined. The two-component plasma is 
known to be equivalent to the sine-Gordon model and 
using this relationship and the results known for this in- 



tegrable field theory, the thermodynamic properties of 
the two-component plasma in the bulk have been exactly 
determined Q in the whole range of stability T < 2. 
However, there are no exact results for confined systems 
for arbitrary Y (with the exception of a Coulomb system 
near an infinite plane conductor p| or ideal dielectric [f| 
wall) . 

It is also well-known for some time that when T = 2 
the sine-Gordon field theory is at its free fermion point. 
This means that the system is equivalent to a free 
fermion field theory and therefore much more informa- 
tion on the system can be obtained. In particular the 
thermodynamic properties and correlation functions can 
be exactly computed even for confined systems in sev- 
eral different geometries and different boundary condi- 
tions 0, H, S E3, O, E3, O O From now on we will 
consider only the case when Y = 2. 

Since at Y = 2 a system of point particles is not stable 
one should start with a regularized model with a cutoff 
distance a which can be the diameter of the hard-core 
particles or the lattice spacing in a lattice model ■ The 
system is worked out in the grand-canonical ensemble 
at given chemical potentials jti+ and /i_ for the posi- 
tive and negative particles respectively. In the limit of 
a continuous model a — > the grand partition function 
and the bulk densities diverge. However the correlation 
functions have a well-defined limit. In this continuous 
limit it is useful to work with the rescaled fugacities 
rn± = 2nde^ fl± /a 2 that have inverse length dimensions. 
The length mT 1 can be shown to be the screening length 
of the system [i| . If external potentials V± (r) act on the 
particles (as in our case) it is useful to define position 
dependent fugacities m±(r) = m± exp(— /?V±(r)). 

Let us briefly review the method of resolution de- 
scribed by Cornu and Jancovici 0, || for the two- 
component plasma. It will be useful to use the complex 
coordinates z = re l6 = x + iy for the position of the par- 
ticles. For a continuous model, a — ► 0, ignoring the possi- 
ble divergences for the time being, it is shown in Ref. 
that the equivalence of the two-component plasma with 
a free fermion theory allows the grand partition function 
to be written as 



2 = det 



2d z 
2d z 



m+ (r) 2d z 
2d z m_(r) 



(2.1) 



Then defining an operator K as 



K 



2d z 



7n+(r) 
m_ (r) 



(2.2) 



2d- z 

the grand partition function H can be expressed as 

5 = dct(l +K) . (2.3) 

The calculation of the grand potential /3fi = — In 3 and 
the pressure p — —dfl/dV where V is the volume (in a 
two-dimensional system this refers to the area), reduces 
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to finding the eigenvalues of K because the grand poten- 
tial can be written as 



n = -fc B T^hi(l + Xi 



(2.4) 



where Xi are the eigenvalues of the operator K. 

On the other hand, the calculation of the one-particle 
densities and correlations reduces to finding a special set 
of Green functions. As usual, the density can be ex- 
pressed as a functional derivative 



p±(r) = m±(r) 



If we define the 2x2 matrix 



(HnS 

6m±(r) 



(2.5) 



as the kernel of the inverse of the operator 

m+(r) 2d z 
2d z m_(r) 



(2.7) 



then the one-body density and two-body Ursell functions 
can be expressed in terms of these Green functions as 

p ai (rx) = m Sl G SlSl (ri,ri) , 
piill(ri,r 2 ) = -m Sl m S2 G SlS2 (r l5 r 2 )G S2Sl (r 2 ,ri) 

(2.8) 

where Si j2 G {+, — } denote the sign of the particles. In 
polar coordinates these Green functions satisfy the fol- 
lowing set of equations 



m+(ri) 
d ri + 



ri - 



a),, 



m_(ri 



G = «5(r 1 -r 2 )I 
(2.9) 

with I being the unit 2x2 matrix. 

The above formalism is very general, it can be applied 
to a variety of situations. In the case we are interested 
in, we will consider two geometries in which the system 
is confined: a disk of radius R and an annulus of inner 
and outer radius i?i and i? 2 respectively. 

The negative particles are supposed to model am- 
phiphathic molecules and therefore they are attracted to 
the boundaries of the system while positive particles are 
not. This is modeled by an attractive external one-body 
potential V- (r) acting on the negative particles near the 
boundary while for the positive particles V+(r) = 0. 

Actually we will consider two models for this potential. 
In the first model (model I) the fugacity m_(r) for the 
negative particles reads, for the disk geometry, 



.(r) = me-P v - (r) =m + aS(r - R) 



(2.10) 



inside the disk, while m+ = m is constant. Outside the 
disk r > i? 2 both fugacities vanish. In the annulus ge- 
ometry, 



inside the annulus. The coefficients a, a\ and a 2 mea- 
sure the strength of the attraction to the walls. In the 
following we will call these coefficients adhesivities. 

In the second model (model II) that we will eventually 
consider the external potential V_ (r) is a step function 
with a range S of attraction near the boundary. This 
model allows us to obtain valuable information regarding 
the frontier regions. Actually we will report here only the 
main results for the annulus geometry with model II in 
section IIV CI further results on this model on the disk 
geometry can be found in Ref. |l5| . 

In the two following sections we will apply the method 
presented here to obtain the grand potential and the den- 
sity profiles of the system. 



III. THE PRESSURE 
A. The grand potential 

To find the pressure of the confined Coulomb system 
we proceed first to compute the grand potential. As 
shown in section ^ the grand potential fl is given by 



/jn = -lnIJ(l + A0 



(3.1) 



where are the eigenvalues of the operator K given 
by equation l|2.2|l . The eigenvalue problem for K with 
eigenvalue A and eigenvector {ip, x) reads 



rn_(r)x(r) = 2Xdgip(r) 
m+(r)V>(r) = 2Xd z x{r) 

These two equations can be combined into 

. . . m+(r)ra_(r) , . 
Ax(r) = +V > x2 X (r) 



(3.2a) 
(3.2b) 



(3.3) 



We now detail the computation of the grand potential 
for the case of a Coulomb system confined in a disk. The 
annulus geometry follows similar calculations. 

We will only use model I for the attracting potential on 
the boundary for the calculation of the grand potential. 
Using model I the fugacities inside the disk read 

m + (r)=m, m_(r) = m + ad(r — R) , (3.4) 

and outside the disk they vanish. Replacing these fugac- 
ities into the eigenvalue problem equations (|3.2|) we find 
that x(r) is a continuous function while ip(r) is discon- 
tinuous at r = R due to the Dirac delta distribution in 
rn_(r). The discontinuity of ip is given by 



i;(R + ,e)~^(R-,8) = jx(R,0)e- w 



(3.5) 



Defining k = m/A, inside the disk r < R, x(r) obeys the 
equation 



rri— (r) 



in 



+ ai S(r - i?i) + a 2 6(r - R 2 ) (2.11) 



Ax(r) = fc 2 x(r) 



(3.6) 
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with solutions of the form 



X (r,9)=A l e M I l (kr) 



and 



iP(r,e)=A l e i( - l -V e I l _ 1 (hr) 



(3.7) 



(3.8) 



where Ii is a modified Bessel function of order I. Out- 
side the disk m_ = m + = and the corresponding so- 
lutions to equations (|3.2[1 are that if) is analytic and x 
anti-analytic, namely, 



X(r) 
V.(r) 



B ie U9 r- 1 
Bie i(l-^ r l-l 



(3.9) 
(3.10) 



In order to have vanishing solutions at r 
essary that 

MR + ,0) = 0, if Z > 1 
X(R, 9) = 0, if I < 



co it is nec- 



(3.11a) 
(3.11b) 



Using these boundary conditions together with the con- 
tinuity of x at R and the discontinuity (|3.5() of ip at R 
gives the eigenvalue equation 



all 



(mR 



/; =0 if i < (3.12a) 

Jj-i = if /> 1 (3.12b) 



The product appearing in equation (|3.1() can be partially 
computed by recognizing that the l.h.s of the eigenvalue 
equation Ij3.12|) for arbitrary A can be written as a Weier- 
strass product 0, 0, Q, 0] . Let us introduce the ana- 
lytic functions 



// +) (z) = [aI l (mzR) + I l - 1 (mzR)](l-l)\ 



(3.13) 
i-i 



mzR 



(3.14) 



By construction the zeros of are the inverse of the 

eigenvalues A for I > and the zeros of f_, ^ are the 
inverse eigenvalues A for I < 0. Furthermore, since 
/ ; (±) (0) = 1, /', (±) (0) = and /j (±) (z) is an even function 
it can be factorized as the Weierstrass product 



where the product runs over all Aj solutions of equa- 
tion 13.12fl for a given I. Then we can conclude that 
the grand potential l|3.1|l is given by 



*(-) 



(-i)+y>// +) (-i) (3.i6) 



i=o i=i 
After shifting by one the index in the second sum and 

rearranging the expression we find the final result for the 

grand potential 



(3.17) 



with 



hw 



(^) 7 "'" /? ' 



(3.18) 



which is the grand potential for a two-component plasma 
in a disk with hard wall boundaries 11 (a = 0) and 



1 + a 



{mR) 



h{mR) 



(3.19) 



is the contribution due to the attractive potential near 
the walls. 

Now we turn our attention to the case of the Coulomb 
system confined in an annulus of inner radius R\ and 
outer radius R 2 . In this case the position-dependent 
fugacity for the negative particles inside the annulus is 
given by 

m_(r) = m + ai$(r - i?i) + a 2 5{r - R 2 ) (3.20) 



The calculation of the grand potential follows similar 
steps as above. One should solve the Laplacian eigen- 
value problem with the appropriate boundary conditions 
given by the continuity of x and the discontinuity of tp at 
Ri and R 2 . After some straightforward calculations the 
final result for the grand potential is 



(3.21) 



// ±) ( js )=n( 1 -^) ( 3 - 15 ) with 



--^ (I^mR^K^mR,) 



I l+1 (mR 1 )K l (mR 2 )) 



(3.22) 
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and 



1=0 
oo 

i=0 



1 + «1 



1 + ^2 



Ii(mR 2 )Ki(mRi) - Ki(mR 2 )Ii(mRi 



I l (mR 2 )K l+1 (mR 1 ) + I l+1 (mR 1 )K l (mR 2 )_ 

Ii +1 (mR 2 )Ki + i(mR 1 ) - Ki +1 (mR 2 )Ii + i{mRi) 
h {mR 2 )K l+1 (mR 1 ) + J m {mR 1 )K t (mR 2 ) 



(3.23) 



The first term f2^ w is the grand potential for a two- 
dimensional two-component plasma confined in an an- 
nulus with hard wall boundaries 1 (a — 0) and the 
second term f2* t is the contribution to the grand poten- 
tial due to the attractive nature of the walls. 

It shoul d be noted that all sums $FTfy . QHty. (|3~22|) 
and (|3.23() above are divergent and should be cutoff to 
obtain finite results. This is due to the fact that the two- 
component plasma of point particles is not stable against 
the collapse of particles of opposite sign for fiq 2 > 2 and 
a short-distance cutoff a should be introduced, a can be 
interpreted as the hard-core diameter of the particles. If 
R is the characteristic size of the system (for instance the 
radius of the disk in the disk case) then l/R is a wave- 
length and the short-distance cutoff a gives an ultraviolet 
cutoff 1/a. Then the cutoff for I (say N) in the sums 
should be chosen of order R/a |illTll|. 



It is interesting to verify if this finite-size expansion 
holds for the systems studied here, in particular if the 
finite-size correction is modified by the special attractive 
nature of the walls considered here. 

Let us first consider the case of the disk geometry. We 
choose to cutoff the sums (I3.18f) and <|3.19ll to a maximum 
value for I equal to R/a and the results given here are for 
a — > 0. The finite-size expansion for large- R of the hard 
wall contribution to the grand potential has already been 
computed in Ref. 0] with the result 

pnZ = -{3p b irR 2 + f3 lhw 2irR+lHmR) + 0(l) (3.25) 

6 



where the bulk pressure pb is given by 

- > 2 



0Pb = -z- 1 + In — 
Zir \ ma 



(3.26) 



B. Finite-size corrections 

It is instructive to study the behavior of the grand 
potential when the system is large. It has been known 
for some time that two-dimensional Coulomb systems in 
their conducting phase have a similar behavior to critical 
systems 0, Il2j . In particular the grand potential of a 
two-dimensional Coulomb system confined in a domain 
of characteristic size L has a large- L expansion 



and the surface tension 7h w for hard walls is given by 

1 1 



/37hw = m 



4 2tt 



(3.27) 



We only need to compute the large-i? expansion 
of H3.19fl . This can be done expressing the Bessel func- 
tion Ii + i(mR) as Ii+i(mR) = I[{mR) + Hi(mR)/(mR), 
using the uniform Debye expansions of the Bessel 
functions valid for large argument 



(3Q = AL 2 + BL 



^lnL 
6 



O(l) 



(3.24) 



similar to the one predicted by Cardy [l7l Il8l | for criti- 
cal systems. The first two terms are respectively the bulk 
grand potential and the surface contribution to the grand 
potential (the surface tension) and are non-universal. 
The logarithmic term is a universal finite-size correction 
to the grand potential, it does not depend on the mi- 
croscopic detail of the system, only on the topology of 
the manifold where the system lives through the Eulcr 
characteristic \- F° r a disk \ = 1 an d for an annulus 

x = o. 



1 Eq. (4.16) of Ref. Illl for the grand potential in an annulus 
with hard wall boundaries is incorrect, however the equation 
above (4.16) is correct and gives our result 13.221 for the grand 
potential 



I,(z) 



27r(/ 2 +Z 2 )V4 
(Z 2 + Z 2)l/4 e „ 



1 



12m z 



3t - 5t 3 
24Z 

9t - 7t 3 
241 



3.28a) 
'3.28b) 



with V = VI 2 + z 2 + I In ( l+ J 2+z2 ) and t = l/Vl 2 + z 2 , 
and using the Euler-McLaurin formula to transform the 
summation into an integral 



N 



TV 



f( x )dx + ±[f(N) + f(0)] (3.29) 



+ ^[f'(N)-f(0)}+--- 

After some calculations, taking first the limit N — > oo, 
keeping only the non- vanishing terms, then taking the 
limit R — > oo and replacing N by R/a, we find that (|3.19|) 
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contributes only to the surface tension giving for the 
grand potential the final result 

[3rP = -PptnR 2 + 2ttR(3 1 + - ln(mi?) + 0(1) (3.30) 
where the surface tension is now given by 



a 171 



a In h 1 — 7r + a 

ma 



I -a 1 



ln(a + 1) 



(3.31) 

We recover as expected the surface tension obtained in 
Ref. for the same system but confined in a slab. 

The universal logarithmic finite-size correction 
(1/6) ln(rai?) is still present and it is not modified by 
the presence of the attractive boundaries. 

For the annulus geometry we are interested in the limit 
Ri — > oo and R 2 — > oo with R2/R1 finite. We proceed as 
above, using also this time the Debye expansion for the 
Bessel functions Gj| 



Kl[z)r 



K[{z)r 



y/2(P + z 2 )V4 



1 - 



3t - 5t 3 



241 



Vg^ + z 2 ) 1/4 e- v 



9t - 7t 3 
241 



(3.32a) 



(3.32b) 



Using l|3.28|l and (|3.32|l one can notice that inside the 
logarithm in l|3.22|l the term Ii + i(mRi)Ki(mR2) is ex- 
ponentially small compared to Ii{mR 2 )Ki + i(mRi) since 
i?2 — R\ — > 00. Also in the contribution from the attrac- 
tive boundaries (|3.23|) the dominant terms are 



N 



N 



1=0 



l + ai 



l + a 2 



KiimRj) 
Ki +1 (mRi) 

Ii + i{mR 2 ) 



Ii(mR 2 ) 



(3.33) 



where i = 1 for the inner boundary and i — 2 for the 
outer boundary. The cutoff N should be chosen ^l| as 
R/a where here R = R 2 x~ x l^~ x > with x = R1/R2 
in order to insure extensivity and recover for the bulk 
pressure pi, the same expression B3.26JI as before. 

In the limit a —> the leading term of the surface 
tension is the same as in the slab and disk geometry 



for a wall with adhesivity a. 

In equation l|3.34|l we do not find any logarithmic 
finite-size correction. There are some terms of the form 
lni?2/-Ri which are order 1 because R2/R1 is finite. This 
is in accordance with the expected formula (|3.24|l for an 
annulus where the Euler characteristic is x = 0. Here 
again the special attractive nature of the walls does not 
modify the universal finite-size correction. 

As a conclusion to this part we might say that the 
logarithmic correction is really universal, not only it does 
not depend on the microscopic constitution of the system 
but also it is insensitive to the existence of a short-range 
one-body potential near the walls. 



C. The disjoining pressure 

Let us first detail the case of the disk geometry. The 
pressure p is given in terms of the grand potential f2 by 



P = 



1 on 

2^RdR 



(3.37) 



Using equation l|3.17|l together with equations (|3.18() 
and H3.19fl gives 



P = Phd + Pat 



(3.38) 



Then, after some calculations we get the final result 
f3n A = -l3 Pb n{Rl ~ R\) + [3 ll 2nR l + /3 72 2^i? 2 + 0(1) 
The surface tension for each boundary is given by 

I -a 2 



m 



, 2N 

cti In h 1 — 7r 

Ri 



(3.34) 

ln(ai + 1) 
(3.35) 



where phd is the pressure for a disk with hard walls 
boundaries given by 



m ^ I l+1 (mR) 
irR Ii (mR) 



(3.39) 



and p a t is the contribution to the pressure due to the 
attractive potential near the walls and it is given by 



/3Pat = 



iR [lf{mR) 



If +1 (mR)] 



{2l + l)I l+1 (mR)Ii(mR) 



2irR 2 ^ 
1=0 



h{mR) [Ii[mR) + aI l+1 (mR)] 



(3.40) 



r 



To study the stability of the system against an external pa = p — Pt defined as the difference between the pres- 
applied pressure one should study the disjoining pressure 
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sure of the system and the pressure of an infinite system 
(the bulk pressure). Let us first consider the case a = 0. 
A proper way to subtract the bulk pressure from expres- 
sion (|3.39|) is by using the equation of state of the infinite 
system j9| 



1 



-n b 



m 

47T 



(3.41) 



where n b is the bulk total density. A simple scaling ar- 
gument shows that for (3q 2 < 2 the equation of state of 
the two-dimensional two-component plasma is poj 



P Pb = [1 - (pq 2 /4)}n b 



(3.42) 



For (3q 2 — 2, the case considered here, the introduction 
of a cutoff a is needed in order to avoid divergences but 
this breaks the scale invariance of the two-dimensional 
logarithmic Coulomb potential giving rise to the anoma- 
lous term (m 2 /47r) in the equation of state. Notice that 
when the cutoff a — > both p b and n b diverge and we 
have Ppd/n b = 1/2 in accordance to the general equation 
of state J2Z21). 

Formally, the bulk density can be written as (see next 
section for details) 



n b 



= — Y Ii(mR)Ki{mR) 

IT ' J 



l=—oo 



m 

7T 



I K +2j2liKi 



(3.43) 



In the above expression and below the omitted argument 
of the Bessel functions is mR unless stated otherwise. 
On the other hand using the Wronskian (2j hKi+i + 
Ii + \Ki = (1/mR) the hard disk pressure (|3.39l) can be 
formally written as 



a m 
PPhd = — 

7T 



i, 

Ll=0 1=1 
Then the disjoining pressure for a — is given by 



(3.44) 



P(Phd - Pb) = /3phd,disj 



m 



E 

.1=0 



hK 



Although the pressure phd and the bulk pressure p b are 
divergent when the cutoff a vanishes, the disjoining pres- 
sure Phd.disj in the case a = proves to be well-defined 
for a — ► and the series ()3.45(l is convergent. 

A plot of the disjoining pressure Phd.disj for a hard wall 
disk as a function of the radius R is shown in Figure ^ 
Notice that Phd.disj is an increasing function of R and it 
is always negative. This shows that in the absence of the 
attractive potential on the walls the systems is always 
unstable for any radius R. This is a common feature of 
the disk geometry with the slab geometry studied in our 
previous work The system without any attractive 
potential on the boundary (a = 0) is naturally unstable. 



<3Phd,di: 



FIG. 1: The disjoining pressure Phd,disj for the disk in the 
case of non-attractive boundaries (a = 0) as a function of the 
radius R. Notice that the disjoining pressure is always nega- 
tive and it is an increasing function of the radius, indicating 
that the system is always unstable. 



For a ^ 0, the disjoining pressure pd is given by 



Pd = Phd,disj + Paid 



(3.46) 



with p a td given by equation (|3.4U|) . Before proceeding 
one should be aware of an important fact. Although the 
first term phd.disj is finite when the cutoff a — * the 
second term on the other hand is divergent when a — > 
0. The sum in equation l|3.4U|l should be cutoff to a 
upper limit N = R/a as it was done in the preceding 
section. Rigorously speaking when a — > the dominant 
term for the disjoining pressure is p a td- For large radius 
R, from last section results (|3.30|) and (|3.36(l on the finite- 
size corrections we can deduce the dominant term of the 
disjoining pressure when a — > 



PPd 



1 

R 



01 



am 2 
=- In 

AttR ma 



(3.47) 



This term is always positive and is a decreasing function 
of R indicating that the system is always stable. 

This is an important difference between the slab geom- 
etry studied in Ref. and the present case of the disk. 
In the slab geometry the disjoining pressure is always fi- 
nite for a — and any value of a. For a slab of width W, 
the large-VF expansion of the grand potential per unit 
area w reads |3j 



- Pb W + 2 7 + O (e- mW ) 



(3.48) 



Then the pressure p = du/dW does not contain any con- 
tribution from the surface tension. On the other hand in 
the disk geometry considered here the existence of the 
curvature makes the surface tension 7 very relevant for 
the disjoining pressure (see equation (|3.47|) ') and since 7 
diverges logarithmically with the cutoff it plays a domi- 
nant role in the stability of the system. 
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FIG. 2: The disjoining pressure pd for the disk as a function 
of the radius R for a = 0.15, 0.21 and 0.3 from bottom to top. 
The cutoff is chosen as ma = 10 -3 . The case a — 0.3 shows 
a stable disk for all radius, in the case a = 0.15 only small 
disks are stable and when a = 0.21 there is the possibility of 
a collapse from large disks to small disks. 



Let us now consider a small but non-zero cutoff a. It 
is expected that the results of our theory in that case 
should be close to the ones of a model of small hard-core 
particles of diameter a. Then there will be a competition 
between the natural unstable behavior of the case with- 
out attractive potential (a = 0) and the attractive part 
to the pressure p a td which is stabilizing. 

Figure [5] shows several plots of the disjoining pressure 
Pd as a function of the radius R for three special values 
of the adhesivity a = 0.15, 0.21 and 0.3 for a fixed cutoff 
ma = 10~ 3 . These three plots show three characteristic 
regimes in which the system can be. 

For large values of the adhesivity, for example the case 
a = 0.3 shown in the Figure|2 the disjoining pressure is 
always positive and a decreasing function of the radius R. 
This is a normal behavior: if the system is compressed 
the internal pressure increases. This indicates a stable 
system for all values of the radius R. 

In the case a = 0.15 (small values of the adhesivity) 
we have two different behaviors of the disjoining pres- 
sure. For radius R smaller than a certain "critical" ra- 
dius i? c , the disjoining pressure is a decreasing function 
of the radius R, indicating again a stable system. How- 
ever, for radius R larger than R c the disjoining pressure 
is an increasing function of the radius R. This behav- 
ior is anomalous, indicating that in this range of radius 
the system is not stable. A very large disk R — > oo is 
marginally stable and will collapse to a disk of smaller 
radius R* where R* is the radius corresponding to pd = 
(see Figure |3J). 

For a = 0.21 (intermediate values of the adhesivity) 
we have a crossover regime. There are now three differ- 
ent behaviors of the pressure characterized by two spe- 
cial radius R^P and R^P . For R < Rc the pressure is 
a decreasing function of R: a stable regime. Then for 
rP < R < R c 2 ^ the pressure is an increasing function 

(2) 

of R: an unstable regime. Finally for R > R c ' the pres- 



FIG. 3: The disjoining pressure pd for the disk as a function of 
the radius R for a = 0.15. The cutoff is chosen as ma — 10~ 3 . 
The radius R c is defined by 5^ = 0. A very large disk R — > oo 
would collapse to a disk of radius R* with disjoining pressure 
Pd{R") = 0. The dashed region for R > R* is not physical. 



. 03 
. 025 

. 02 
. 015 

0.01 

!t(2) 

. 005 



R<J> 



FIG. 4: The disjoining pressure pd as a function of the radius 
R for a = 0.21 and a cutoff ma = 10~ . The theoretical result 
for the pressure shows a non-physical region (dashed line). 
The correct pd vs. R curve can be obtained by a Maxwell 
construction, thus showing that in this case there is a first 
order transition (collapse) from a large disk (radius Rb) to a 
small disk (radius R a ). 



sure is again a decreasing function of R: another stable 
regime. One can do a Maxwell construction do determine 
the correct dependency of the pressure against the radius 
as shown in Figure 0] In this regime there is a first order 
transition (collapse) from a large (but finite) disk (radius 
Rb) to smaller disk with radius R a (see Figure H}. 

Two of the three regimes illustrated here, small adhe- 
sivity (for example a = 0.15) and large adhesivity (for 
example a = 0.3), also occur for the slab geometry stud- 
ied in Ref. |3j. In that case they were separated by the 
special value a = a c = 1. The disk case considered here 
is however more rich since there is also a crossover regime 
between the two, for intermediate values of a (for exam- 
ple a = 0.21), with the possibility of stable large radius R 
disks, a forbidden (unstable) intermediate range of radius 
R and then again stable small disks. 

The values of a characterizing the different regimes 
are however highly dependent on the cutoff and so are 
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the values of the "critical" radius i? c , Ri' and .Re- 
values a — 0.15, 0.21 and 0.3 of Figure |21 are for a cutoff 
ma = 10~ 3 , for other values of the cutoff these values 
will change. For this reason the analysis above should be 
considered with some care. For very small values of the 
cutoff a we expect that the results obtained here should 
reproduce asymptotically those of a primitive model of 
hard-core particles with radius a. 

To conclude this section let us mention that in the 
annulus geometry the situation will be similar to the disk 
one. In that case one can consider the pressure on the 
inner boundary or the pressure on the outer boundary. 
For both pressures the dominant term when the cutoff 
a — ► will be given in terms of the surface tension 7 
as in equation (|3.47l) . Then in the strict limit a — > 
the annulus geometry will be stable if a > as in the 
disk case. Considering a small but non-zero cutoff a will 
lead to a similar discussion as in the disk with different 
regimes some of them exhibiting a collapse. 



The Equation (|4.2|l has solutions of the form 



G±±(r 1; r 2 ) = £ J2 e* m - e2) [Ilimr^Kiimr^ 



l=-c 



Affa)Ii(mri)] 



(4.4) 

max(ri,r2). The re- 



where r< = min^i,^) and 
maining Green functions can be found using Eq. (|4.3|) . 

If ri is outside the film while r 2 is fixed inside the film 
m+(ri) = m_(ri) 
are 



G +T (ri,r 2 ) 
G_ T (n,r 2 ) 



0, the G SlS2 that satisfy Eq. (HI 
GfaMne" 1 ) 1 , (4.5) 

I — — OO 
+ OO 

J2 D l (r 2 ,9 2 )(r 1 e-^)~ l . (4.6) 



IV. THE DENSITY PROFILES 

In this section, we study the densities of the Coulomb 
system confined inside a disk and inside an annulus. The 
densities can be obtained by computing the Green func- 
tions introduced in Sec.[n] First we will solve and exten- 
sively study model I for the disk then the annulus. Then 
we will briefly consider the results that can be obtained 
using model II. 



So, in order to have finite solutions at r± = 00 it is nec- 
essary that for n > R 



Ci(r 2 ,e 2 ) = for^>0 
A(r 2 ,0 2 ) = O for^<0, 



(4.7) 



Eqs. (|4.1|) and 14. 7|) are the boundary conditions that 
complement the differential equations (|4.2[1 and (|4.3|) for 
the Green functions. 

Solving for the coefficients Af , we arrive at the follow- 
ing expressions for the Green functions for < ri, 2 < R: 



A. The disk 

In the disk geometry, for model I, the fugacity m_ (r) is 
given by Eq. i|2.1U[l . From Eq. I|2.9|) . when r 2 ^ R, we see 

that G__ and G |_ are continuous for n = R. However 

because of the Dirac delta distribution in the definition of 

m_(r) the functions and G++ are discontinuous at 

r\ = R. The discontinuity can be obtained from Eq. (|2.9|) 
if r 2 ^ n: 

G +T ( ri =R-,r 2 ) - G +T {r 1 ^R+,r 2 ) 

= a e- i9l G_ T (n = R,r 2 ) 

(4.1) 

If both points ri and r 2 are inside the disk but not on 
the boundary then Eq. 1)2. 9fl lead to a Hclmoltz equation 
for G++ and for G__ 

[m 2 - A] G±±(r 1 ,r 2 )=m5(r 1 -r 2 ) (4.2) 

and the other Green functions can be obtained from 

e ±idl ( i \ 

-dry T — d Sl G±± (ri,r 2 ) = G T ± (ri,r 2 ) 

m \ r\ J 

(4.3) 



G++(ri,r 2 ) = — K (m\ri - r 2 | 



m 
2^ 



1=0 

aK l+1 - Ki 



Ki 



Ii + i(mri)Ii + i(mr 2 ) 



h 



Ii(mri)Ii(mr 2 ) 



(4.8) 



and 



G — (ri,r 2 ) = — K (m\ri -r 2 | 

Z7T 



m 
2^ 



1=0 



Ii(mn)Ii(mr 2 ) 



Ki 



oih+i + h 



Ii+i{mri)Ii + i(mr 2 ) 



(4.9) 



with // meaning // (mi?) , and the same convention for the 
Bessel function Ki. 

The one-particle densities are given in terms of these 
Green functions as 



p B (r) = m s (r)G ss (r,r) 



(4.10) 



As we explained in section^ the continuous limit model 
presents divergences in the expressions of the densities. 
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This is seen in the term K^rnr 12) that diverges logarith- 
mically as r 12 = |i"i — r 2 | — * 0. So we impose a short 
distance cutoff a. One can think that particles are disks 
of diameter a, so the minimal distance between particles 
is a. 

The first term in Eqs. I|4.8|) and (|4.9jl gives the bulk 
density pb, the density of the unbounded system as cal- 
culated in Ref . [!j ■ For a — > 



Pi 



9 9 

m . . to 
Pb = Pb = —K (ma) ~ — 



In- 



'7 



(4.11) 



where 7 ~ 0.5772 is the Euler constant. In the second 
terms of Eqs. I|4.8[) and (|4.9|l the sum can eventually di- 



verge when ri = r 2 = r for certain values of r (in the 
boundaries) so we should impose a cutoff \l\ < N = R/a 
as it has been done in the expressions of the pressure and 
grand-potential obtained in the last section. 

Because of the form (|2.10|l of m_(r), the negative den- 
sity can be written as 

p-{r) = (l + -5(r- R))p*_(r) (4.12) 
\ m J 



where p*_ can be seen as the density of non-adsorbed 
particles. For the positive particles p + — p* + . Finally we 
have 



P%(r) 
P-(r) 



Pb 



Pb 



rn 
2 



2 00 



aK l+1 - Ki 



"2T ^ 



Ki o, , 
— 7 Z (mr) + 



aKi+i - Ki j2 
ah+i + h 



(4.13a) 
(4.13b) 



The non-adsorbed charge density p* = p* + — p*_ can be 
obtained from the above expression and using the Wron- 
skian of the Bessel functions J;A'; +1 + Ii +1 K[ = 1/mR, 

P \T)=^f (4.i4) 

P{ > 2ttR^ (ali+i + li)li 1 ' 

Finally, the total charge density 

P( r ) = P+(r) - P-(r) 

= P+(r)-(l + -S(r-R))p*_(r) (4.15) 

V TO / 

= p*{r) -aS(r-R) 

has a non-adsorbed part p* (r) and a adsorbed "surface" 
charge density in the boundary 



<7_ =-p*-{R) 
m 



This surface charge density cr_ comes from the S(r — R) 
part of the negative charge density. Writing formally 
the bulk part of the density as pb = (to 2 /27t) ^2i e ^hKi 
one can obtain the following expression for er_ from 
Eqs. gT^I and ^4.13bt 



\ ^ 



2vri? f-^ a/i+i + /; 



(4.17) 



Actually the adsorbed charge density should obey two 
special relations. The first is a sum rule that expresses 
the global electro-neutrality of the system 



Ra- 



p* (r) r dr 



(4.18) 



(4.16) Using the indefinite integral 



rn R 



x(lf +1 (x) + lf(x))dx 



{mRf 



[lf +l - hh+2 + If - h-ih+i] = mRIiIi+i 



(4.19) 



obtained using the recurrence relations [2l| for the Bessel 
functions J; , the sum rule l|4.18(l is immediately shown to 
be satisfied. 

On the other hand the adhesivity a can be thought as 



a sort of fugacity that controls the number of adsorbed 
particles and one can obtain the total number of ad- 
sorbed particles 2nRa- from the grand potential £1 by 
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p(r) /m z 
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the ones that lie in the boundary (the surface charge cr_). 
Therefore, the negative density p-(r) decreases near the 
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FIG. 5: The charge density profile p(r) for a disk of radius 
R — 10/m. From bottom to top the adhesivity a — 0.25, 0.5, 
0.75. 



using the usual thermodynamic relation 



2ttR<7- 



-a f3 



da 



(4.20) 



This relationship is also immediately shown to be satis- 
fied from the expression (|3.19l) for fl^ t , the part of the 
grand potential that depends on a. 

For a large disk, R — > oo, the dominant part of the 
grand potential that depends on a is the surface tension 
7 given by equations l|3.31|l and (|3.36() . Then we have 



-/3 a 



d-j 

da 



(4.21) 



a relation already shown to be true in Ref. 3] for the 
same system near a plane attractive hard wall. Using 
Eq. i|33T)l into Eq. (|QT|) gives explicitly 



-in 



a In ■ 



1 



ma 



ln(a + l) + l 



(4.22) 



Thus recovering a known result from Ref. 0. Notice 
that for large disks the adsorbed surface charge density 
becomes independent of the radius R. This result can of 
course by also obtained directly from Eq. (|4.17l) in the 
limit of large- R using the Debye expansions (|3.28|) of the 
Bessel functions. 

In Figure El we can see the charge density profile p(r) 
(actually only the non-adsorbed part p*{r) is shown). It 
can be seen that the negative surface charge — cr_ (not 
shown in the figure) is screened by a positive layer near 
r = R. As the adhesivity is increased, cr_ increases, 
and therefore the density of the screening positive layer 
shown in Figure increases. This screening layer has a 
thickness of order of magnitude to" 1 as expected since 
to" 1 is a measure of the screening length [9j. 

In Figure|H| we can see the distribution of the negative 
particles (two lower curves) and positive particles (two 
upper curves). The negative particles feel repelled by 



FIG. 6: The density profiles p+(r) and p-(r) for a disk of 
radius R — 10/m. The bulk density has been subtracted from 
both densities. The upper curves correspond to the density 
of the positive particles p+(r) and the bottom curves to the 
density of the negative particles p-(r). The adhesivity a — 
1 for the topmost curve (p_|_(r)) and the bottommost curve 
(p_(r)). The two other curves correspond to a = 0.5. 



frontier. The positive charge density p+{r) inside the 
disk can be observed in the upper curves of Figure El 
Due to the large accumulation of negative charges at the 
frontier, the positive particles tend to shield this charge. 
Hence, there is an increase of the positive density near 
R. 

As we increase the adhesivity, the negative adsorbed 
charge on the border gets larger while on the inner region 
there is a stronger repulsion of the negative charges and 
an increase in the number of positive particles shielding 
the outer charge. 



B. The Annulus 

For the annulus geometry with a delta distribution 
modeling the external attractive potential (model I), we 
follow the same reasonings as for the disk geometry. The 
solution to Eq. (|4.2|l for R\ < r < R 2 is of the form 

G±±(ri,r 2 ) = £ J2 elWl2 ft(™-<)Jr,(ror>) 

I— — OO 

+ Af(r 2 )Ii(mri) + B^(r 2 )Ki(mri)] 

(4.23) 



2. As in the disk case, the 
. are continuous while G++ 



where we defined #12 = 61 - 

Green functions G and G 

and G-j are discontinuous at r = Ri and r = R 2 with a 

discontinuity given by Eq. (|4.1(l replacing the adhesivity 
a by Qfi for the discontinuity at r = Ri and by a 2 for 
the discontinuity at r — R 2 . 

Solving for the coefficients we finally find 
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vn 

G__(n,r 2 ) = — KbMn - r 2 |) (4.24a) 
- ^E n (i)' Ii(mn)Ii(mr 2 ) (4.24b) 

-«0 12 d(1) 

" 2^E n (i)' ^(mn)^(mr 2 ) (4.24c) 
to °° e~^ 12 C^ 

+ ttE" — ^[/iHj^HJ + ^K)/,^)] (4.24d) 
2 ^ «=Z D i 

m ~ 6^ + 1)612^(2) 

+ ^E ^/i+iK)^!^) (4.24c) 

°° e i(i+i)ei2 r(2) 

+ -^-^(mnj^+^mra) (4.24f) 

i=o D t 
m ~ e *C+i)0i2c( 2 ) 

~ 2~E (2) ' [^+i(mri)7j + i(mr 2 ) +/; + i(mr 1 )Jf l+1 (mr 2 )] 

^ Z=o A 

(4.24g) 



and 



G ++ (ri,r 2 ) = — tf (m|ri - r 2 |) (4.25a) 

Z7T 

+ ^E n( i) ' Ii+i(mn)I l+1 (mr 2 ) (4.25b) 

m ~ e -»('+i)ei25( 1 ) 

+ ^E n (i) ' Jft + i(mri)ifr + i(mr2) (4.25c) 

+ 7T / , m — [7i+i(mri)/<( + i(mr 2 ) + ^i + i(mri)/| +1 (mr 2 )] 

(4.25d) 

°° iiei 2 4(2) 

- ^E n(2 )' Ji(r»ri)J,(mr 2 ) (4.25e) 

oo U0 12 d(2) 

- f E 6 (2/ ^(mrO^Cmra) (4.25f) 

oc giJ0 12 ^»(2) 

- £E~ ^— [^(^1)^(^2)+ #i(«wi)Ji(mr 2 )] (4.25g) 
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with 



.4 



B 



(i) 



D 



(i) 



.4 



(2) 



B 



(2) 



(2) 



D 



(2) 



= K 



r(2) 



(1) 



(1) 



K. 



(2) 



= I, 



(2) 
i 



/v. 



(1) 

(1) 



ail, 
(a 3 K 



I- K, 

'i+i 
7 (i) 



(i) 
i+i 



(i) 



A', 



(i) 



(2) 
i+1 



= I 



1 + 1 

(1) 
/ + 



r(2) 



r (2) 



(2) 



(1) 
2+1 



a 2 K l 



(2) 
2+1 

(2) 
2+1 



-K, 



(2) 



A" 



(2) 



(2) 



A' 



(i) 
2+1 



(2) 



(4.26a) 
(4.26b) 
(4.26c) 
(4.26d) 
(4.26e) 
(4.26f) 
(4.26g) 
(4.26h) 



where 7 ; = Ii(mR\), l\ — Ii(mR 2 ) and the same 
convention for the other Bessel functions. 

The individual densities are obtained putting ri = 
r 2 = r in the above expressions since p s (r) — 
m s (r)G ss (r,r). The above expressions are quite long, 

however one can notice some symmetries between G 

and C++. Beside a change of sign and the phase factor 
(which is irrelevant in the calculation of the densities) 

the coefficient of the term Ii(mri)Ii(mr 2 ) in G is the 

same as the one for Ii + i(mr{)Ii+i(inr 2 ) in C++ and so 



on. The charge density can be written again as 

p(r) = p*(r) - cr a) <5(r - R x ) - a { _ ] S(r - R 2 ) (4.27) 
with a non-aborbed part p*(r) and a surface density 

<y { - 2) = (a 1 , 2 /m)p*_{R 1 , 2 ) (4.28) 

of adsorbed negative particles in i?i and R 2 respectively. 
The non-adsorbed charge density is 



2 00 



2=0 



(Cft + Cl\) [l? +1 (mr) + If(mr)} + (C«f + C« K ) [K? +1 (mr) + K?(mr)} 



+2(Cif + Cif)[Ii+i(mr)Ki +1 (mr) - Ii(mr)Ki(mr)} 



(4.29) 



with the coefficients that depend on ot\ given by 



Either way the result is the same as expected 



(i) 



D 



(i) 



KK 



B 



(1) 



D 



(i) 



IK 



c, 



(1) 



D 



(i) 



(4.30) 



and the ones that depend on a 2 are 



C. 



ii 

2.1 



n KK _ 



D 



(2)' ^2,i 



B 



(2) 



n IK _ 
J °2.( — 



D 



(2) 



(2) 



Dp) 



(4.31a) 



The adsorbed charge density in each boundary can be 

computed by replacing G (Ri,2, Ri,2) from Eq. (|4.24() 

into Eq. or by using the thermodynamic relation 



(2) 



2ttAi 



E 

1=0 



If'Kf' - Kf'I, 



(2)r(l) 



D 



(i) 



(4.33a) 



oo j(2) K (l) _ K (2) r (l) 

« 2 ^VVl_|^ (4 33b) 



2ttA 



1=0 



D 



(4.32) Using the indefinite integrals (|4.19jl and |2S 
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niR 



■m R 

(K?(x)+K? +1 (x))xdx 

(Il+i(x)Ki + i(x) - h{x)Ki(x))xdx 

I 



-mRKi+iKi 



(4.34) 
(4.35) 



p (r) /m 
. 3 
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FIG. 7: The charge density profile p(r) for an disk of inner 
radius Ri — 1/m and outer radius R2 = 5/m. The adhesivi- 
ties in each boundary have been chosen equals ai = Q2 = a. 
The upper curve correspond to a = 1 while the lower one to 
a = 0.25. 



FIG. 8: The difference between the surface charge in the inner 
boundary and the surface charge in the outer one — u^ 2 ' 
as a function of a for an annulus of inner radius Ri — 1/m 
and outer radius R2 = 5/m. 



one can verify that the electroneutrality sum rule 



fR-2 

/ p* (r) rdr = R 1 a W + R 2 <J { ^ 

JRx 



(4.36) 



is satisfied. For the reader interested in the details of 
this calculation, an extented version of this manuscript 
is available online |23j|. 

Figure shows a plot of the charge density p(r) for 
an annulus with inner radius R\ — 1/m and outer ra- 
dius R 2 = 5/m for ai — a 2 — 0.5 and ai = a 2 = 1. 
In both figures on can see a positive layer of charge 

in each boundary screening the adsorbed negative sur- 
(1 2) 

face charge densities a_ . As a\ t2 increases the ad- 
sorbed charge increases and so does the positive layer. 
Although it is not perfectly clear in Figure there is 
actually slightly more adsorbed surface charge density 
in the inner boundary that in the outer. This can be 
seen in Figure that shows the difference between the 
adsorbed surface charge in the inner boundary and the 



surface charge in the outer one, <r_ — &— J , as a function 



r (2) 



of a. Figure [S] clearly shows that a\ = a 2 . 

On the other hand the total charge on the inner boundary 
is smaller that the total charge on the outer boundary: 
2irRia_ < 1-kR 2 ol\ This can be seen directly from 
Eqs. (|4.33|) . If 0:1=0:2 = a, each term in the series of 



the difference R\ erf 1 ' — R 2 from Eqs. I|4.33|l is of the 



(2) 



form 



o 



(b, 



J l+D 



(4.37) 



where and are the denominators in each term 
of the sums in Eqs. (|4.33a|l and Ij4.33bjl respectively, and 



(2)r(l) 



a-i, 



The sequence (6;)/gn has the 



property of being monotonically increasing with I. Then 
we conclude i?i < R 2 



C. Model II 

We now briefly consider model II and some of its results 
in the annulus geometry. 

In the annulus geometry, there are two outer regions, 
each one has a width 5, and an inner region that has a 
thickness R 2 — Ri — 25. We can distinguish three regions 
for the annulus geometry: the inner border Ri < r\ < 
Ri+S (region 1 ) , the bulk of the film Ri+S < r% < R 2 — 5 
(region 2) and the outer border R 2 — 5 < n < R 2 (region 
3). So for this model in the annulus geometry, we define 
the position dependent fugacities as 



m+(n) 
m_ (n) 



m2 if r\ G region 1 
m if 7*1 G region 2 
m 2 if n € region 3 



(4.38a) 
(4.38b) 
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The fugacity in the border regions is m 2 = 
mexp(— (3U-) where C/_ < is the value of the external 
potential V- (r) near the boundary. Notice that m 2 > m. 
It is clear that model I is the limit of model II when 5 — > 
and m 2 — > oo with the product to 2 (5 = a finite. 

For this model II, the following symmetrization of the 
problem explained below will be useful. If an external po- 
tential is acting differently on positive and negative par- 
ticles, as in our case, in order to symmetrize the problem 
for the two types of particles, Cornu and Jancovici 
propose to define m (r) and V (r) by the relations 



p (r) /rrr 



1 



. 5 

■0.5 
-1 



1 1 
1 \ 










3 


4 





m s (r) = to (r) exp [-2sV (r)] , s = ±1 (4.39) 

where V(r) in general is different from the external V±(r) 
potentials. Knowing that V_(r) is constant in each re- 
gion, in terms of the modified Green functions 

5siS2 (r 1 ,r 2 )=e- SlV ( ri )G SlS2 (r 1 ,r 2 )e- s ^( r2 ) (4.40) 

the system of equations (|2.9() can be shown to be equiv- 
alent to 



(m(ri))' 



—T-r -d ri T -de, ) g±± (ri,r 2 ) = g T ± (ri,r 2 ) 
m(ri) \ n 



where A is the Laplacian operator and 



m(r) 



too = (toto 2 ) 1 / 2 if r in regions 1 or 3 , 
to if r in region 2 . 

(4.43) 

and the potential defined in Eq. 1)4. 39|) is given by 



exp (V(r)) 



(to 2 /to) 1 ^ 4 if r in regions 1 or 3 , 
1 if r in region 2 . 

(4.44) 

To obtain the Green functions g±± we solve Eq. (|4.41|) 
which is a standard inhomogeneous Helmoltz equation in 
each region: inner border R\ < r < R\ + 5 (region 1), 
bulk R\ + 5 < r < i? 2 — 5 (region 2) and outer border 
R—6 < r < R (region3). Outside the annulus r < Ri and 
r > i? 2 the fugacities vanish and the equation satisfied by 
G±± is the Laplace equation. The Green functions G SlS2 
are continuous at each boundary r = i£i, r = R% + 6, 
r = i? 2 — 6 and r — i? 2 . Remember that here the Green 
functions G SlS2 are different from the auxiliary functions 
g Sl s 2 - Their relationship is given by Eq. (|4.4U|I . While 
the Green functions G SlS2 are continuous, the auxiliary 
Green functions g Sl s 2 are m general discontinuous at the 
interfaces. 

In the regions of interest, the form of the solution is, 



FIG. 9: The charge density profile for an annulus with in- 
ner radius Ri = 1/m, outer radius R2 = 5/m in model II. 
The width of the borders are S = 0.5/m. The fugacity in the 
border regions is: for the solid line curve rri2 = 16m corre- 
sponding to a value of mo = 4m and for the dashed line curve 
m,2 = 4m corresponding to a value of mo = 2m. The cutoff in 
the sums of Bessel functions in the expression of the density 
has been chosen as N — 100. 



in general, 



0±±(ri,r 2 ) =m(ri)<S(ri -r 2 ) (4.41) S±±(ri,r 2 ) 



m(ri) 
2tt 



(4.45) 



J2 e u ^ s ^ [/,(m(n)r<)lf,(»n(»-i)r>) 

! — — 00 

+ S i ± (r 2 )/Km(r 1 )r 1 ) + C i ± (r 2 )K i (TO(r 1 )r 1 )] 



(4.42) with m(r) = mo in the border regions 1 and 3 and 
m(r) = to in the bulk region 2. Imposing the conti- 
nuity of the Green functions G SlS2 gives a linear system 
of equations for the coefficients B-f and Cf 1 which can be 
conveniently solved with the aid of matrix algebra soft- 
ware like Mathematica. 

The explicit expression of the coefficients obtained 
from the solution of the linear systems ^| are to long 
to reproduce here, however we show in Figure [5] a plot 
of charge density profile for two values of to 2 , to 2 = 16m 
and m 2 = 4m. 

The discontinuity of the charge density profile near 
the interfaces is due to the fact that the density pro- 
file of the negative particles is discontinuous there. This 
is expected for the negative particles since the poten- 
tial V-(r) is discontinuous across the interfaces. It is 
clear from our formalism that the density profile should 
have the same kind of discontinuity that the Boltzmann 

factor of V-(r) since p~(r) = m_(r)G (r, r) with 

m_(r) = mexp(— /3V_(r)) and G continuous at the 

interfaces. Actually this is a more general result of statis- 
tical mechanics of inhomogeneous fluids at discontinuous 
interfaces. In general for a fluid with density n(x) near 
a planar interface characterized by an external potential 
Vext(x) eventually discontinuous at x = the y- function 
exp(pV c ^t(x)) n(x) is continuous pij . 

There is a higher density of particles (both negative 
and positive) in the borders that in the bulk: for higher 
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values of the fugacity m% in the border (the external po- 
tential V_(r) is more attractive), the density in the bor- 
ders is higher. In region 2, the density far from the in- 
terfaces is close to the bulk value (I4.11|) . In the borders, 
away from the interfaces the density tries to be close to 
the new bulk value given by Eq. I|4.11|) replacing m by 
mo- If the width S of the border region was much larger 
than tUq 1 , one would expect that the density far away 
from the interfaces converges to a bulk value according 
to the value of m(r) in that region. This is true for both 
the negative and positive of particles. Due to the natu- 
ral tendency of the system to be electrically neutral, the 
positive particles try to follow the negative ones so the 
system is not locally charged. However at the interfaces 
r = Ri + 5 and r = R2 — 5 there remains a non-neutral 
charge density that can be seen in Figure Actually 
one can see in Figure ||a double charged layer. Inside 
the border region (1 or 3) there is a negative charge den- 
sity layer and outside the border, in the bulk region 2, 
a positive layer, the same one that was previously ob- 
served with model I in Figure These double layers are 
concentrated near the interfaces at r = Ri + S and at 
r = i?2 — 5 and have a thickness of order mT 1 for the 
positive layer in region 2 and m^ 1 for the negative layer 
in regions 1 and 3. 

V. SUMMARY AND PERSPECTIVES 

The present solvable model studied here gave us in- 
teresting information about the behavior of confined 
Coulomb systems with attractive boundaries. This sys- 
tem has an induced internal charge on the boundary 
which is created by an external potential which is not 
of electrical nature. This potential only acts on the neg- 
ative particles, while the positive particles are unaffected. 
This is not the usual situation that has been studied ex- 
tensively in the past, where the system is submitted to 



electrical forces due to possible external charges. 

First, we found that large systems exhibit the same 
finite-size corrections that for systems without attrac- 
tive boundaries, confirming again the universal nature 
of these finite-size corrections. Studying the disjoining 
pressure we found that the attractive boundaries have a 
stabilizing effect. This was noticed also in our previous 
work Q, however the curvature in the present case is 
very important. It makes the surface tension to be the 
predominant contribution to the disjoining pressure, as 
opposed to the slab geometry. Then, we conclude that 
the curvature has also a stabilizing effect on the system 
in comparison to the slab geometry in which the system 
can be unstable for low values of the adhesivity. 

The study of the density profiles gives information 
about the structure of the system. As expected, there 
are some adsorbed charges on the boundary and these are 
screened by a positive layer of charge inside the system. 
We were able to check explicitly an electro-neutrality sum 
rule and a few relations that the adsorbed charge in the 
boundary satisfy. 

It would be interesting to know what features of the 
present model are universal and which are not. A step 
toward answering this question can be obtained by study- 
ing another solvable model of Coulomb system, the one- 
component plasma. A preliminary study of this system 
was done in Ref. 0] and this will be the subject of a 
future paper. 
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